Velocity-correlation distributions in granular systems 
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We investigate the velocity-correlation distributions after n collisions of a tagged particle un- 
dergoing binary collisions. Analytical expressions are obtained in any dimension for the velocity- 
correlation distribution after the first-collision as well as for the velocity-correlation function after 
an infinite number of collisions, in the limit of Gaussian velocity distributions. It appears that 
the decay of the first-collision velocity-correlation distribution for negative argument is exponential 
in any dimension with a coefficient that depends on the mass and on the coefficient of restitu- 
tion. We also obtained the velocity-correlation distribution when the velocity distributions are not 
Gaussian: by inserting Sonine corrections of the velocity distributions, we derive the corrections to 
the velocity-correlation distribution which agree perfectly with a DSMC (Direct Simulation Monte 
Carlo) simulation. We emphasize that these new quantities can be easily obtained in simulations 
and likely in experiments: they could be an efficient probe of the local environment and of the degree 
of inelasticity of the collisions. 

PACS numbers: 05.20.Dd, 51.10.+y, 02.50.-r 



I. INTRODUCTION 

The dynamics of hard-core particles consists of suc- 
cessive binary collisions. For atomic systems, the equi- 
librium state can be reached, and is characterized by 
velocity distributions which are purely Gaussian. Con- 
versely, in the presence of dissipation, i.e. for granular 
particles, no equilibrium exists, but when an external 
source of energy is present the system reaches a steady 
state whose properties can be compared to the equilib- 
rium state of atomic systems [lj]. At low to intermediate 
densities, spatial correlations are not responsible for non- 
Gaussian deviations [H- 

The short-time dynamics is usually analyzed by means 
of the velocity autocorrelation function. This quantity 
provides an average of the scalar product of the velocity 
at time with the velocity at time t. The characteris- 
tic time of this correlation function corresponds to the 
time needed for the system to lose memory of the ini- 
tial configuration of the velocities. Some progress has 
been made recently by investigating the collision statis- 
tics: Visco et. aZ0, [j showed that the free flight time 
distribution is not exponential, even in the low-density 
limit (such a behavior was observed in Molecular Sim- 
ulation of hard spheres some years agoQ). Deviations 
from the Poisson law of the number of collisions can be 
captured in the framework of the Boltzmann equation 
and agree with molecular simulation results. 

We introduce here a new quantity by focusing on colli- 
sion events irrespective of the time when individual col- 
lisions occur. We consider the scalar product between 
the velocity before a given collision and the velocity af- 
ter n collisions. (Note that this quantity is distinct from 
the distribution of velocity for a hard sphere on collision, 
which characterizes the distribution of the relative veloc- 
ity of colliding spheres, for which Lue[|| obtained exact 
results in three, four and five dimensions). The informa- 
tion obtained is not only the average of the scalar prod- 



uct, but the full distribution of the scalar product be- 
tween the velocity before and after a sequence of n colli- 
sions. When n = 1, this corresponds to the probability of 
the scalar product between the pre- and post-collisional 
velocities during a collision; it is worth noting that the 
first moment of the distribution does not correspond to 
the velocity correlation function at the mean collision 
time: indeed, the probability distribution is built for col- 
lisions occurring at different collision times, whereas the 
velocity correlation function corresponds to the velocities 
scalar product at a given time. 

As the number of collisions increases, the correlations 
between velocities decrease and the distribution evolves 
progressively to the asymptotic form where the velocities 
are uncorrelated. The paper is organized as follows: in 
section II, we obtain the probability distributions in the 
limit of an infinite number of collisions, i.e. when the 
velocities are completely uncorrelated, in any spatial di- 
mension. In section III, we derive the first-collision veloc- 
ity distribution in any dimension. Section IV is devoted 
to the corrections induced by the non-Gaussian behav- 
ior of granular systems, and we compare the analytical 
results to DSMC results. Velocity-correlation distribu- 
tions at the second and higher collisions are analyzed in 
Section V. 

Let us define, right away, the central quantity of this 
study, the velocity-correlation distribution P n (z) at the 
nth collision : 



P n (z) = (S(z 



;)) 



(i) 



where the brackets denote a statistical average in a given 
steady state, v denotes the precollisional velocity of a 
tagged particle before the first collision and v* the post- 
collisional velocity of the same particle after the nth colli- 
sion. Wc only consider the case of homogeneous systems. 
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II. VELOCITY-CORRELATION 
DISTRIBUTION IN THE INFINITE COLLISION 
LIMIT 

We first consider the situation where the number of 
collisions is very large, i.e. the velocity before the first 
collision and the velocity after a large number of collisions 
become uncorrelated. The probability distribution of the 
scalar product Pao{z = v.v*^) is then given by 



PM = 



dvdv^ ) /(v)/(v* 00 )^-v.v* 



(2) 



Let us consider the generating function P oc> (k) 
J dke lkz Poo(z). One has 



P^(k) 



dvdv* 00 /(v)/(v* 



(3) 



We assume that the velocity distribution can be factor- 
ized as /(v) = n«=i f( v a) where a is an index running 
over all components of the velocity and d the space di- 
mension. The Cartesian components of the velocity are 
independent random variables and the generating func- 
tion is then also the product of the generating functions 
of each Cartesian component: 



P£Xk) = {P£Xk)) d , 



(4) 



where P^\k) is the generating function for the one- 
dimensional problem. 

Poo (z) can be obtained from Eq. ([2]), which gives 



dv 



It is worth noting that PoW is the Mellin convolution 
of the two velocity distributions as this distribution is 
that of the product of two independent random variables. 
When the velocity distribution f(v) is Gaussian 



f(z/v)f(v). 



(5) 



f(v) 



M 



(6) 



where M and T are the mass and the temperature of the 
tagged particle, respectively. 

can be explicitly obtained and is then equal to 



ttT 



(7) 



where K$(z) is the modified Bcssel function of second 
kind. 

The basic property of P^ is that the distribution is 
symmetric because the velocities are uncorrelated. The 
behavior of P^ is intriguing at small values of z, as one 
observes a logarithmic divergence at z = 0. This means 
that there is an overpopulatcd density of very small scalar 
products even though the velocity distribution remains 



finite for very small velocities. For large velocities, Pib 
decays as exp(— \z\M/T)/y/z for large values of z, i.e. 
less rapidly than the original velocity distribution which 
has a Gaussian decay. 

In two and more dimensions, the velocity-correlation 
distribution can be obtained by noting that the Fourier 
transform of Eq.([7]) (or integrating Eq.Q), leads to 



i 



(8) 



By inserting Eq.© in Eq.Q, the generating function 
Poo\k) in d dimensions is then 

pW(fc) = (i + fc 2 r d/2 - (9) 

The inverse Fourier transform can be calculated in any 
dimension: in 2D. the Fourier transform has a Lorentz 
profile, which gives in real space 



p£Hz) = 

and in three dimensions 

P£H*): 



Mexp(-\z\M/T) 
2T : 



(10) 



(11) 



In 2D and in 3D, the probability distribution is no longer 
singular at the origin. However, there exists a non- 
analytic behavior which is a cusp in 2D and a cusp in 
the derivative in 3D. 

For completeness, the solution in odd dimensions; 
Pco\z) is given by 



(d-2)! 



2(d-3)/2(( d _3)/ 2 )! 



K((d-l)/2,\z\) 
(12) 



collision where K((d — l)/2, |z|) is the modified Besscl 
function of second kind of order (d— l)/2. 



III. FIRST-COLLISION 
VELOCITY-CORRELATION DISTRIBUTION 

In order to have tractable expressions for the first- 
collision velocity-correlation distribution Pi(z), we as- 
sume that the the "molecular chaos" is valid, i.e. that 
there are no correlations between the particles before col- 
lision. The joint velocity distribution of the tagged par- 
ticle and the bath particles is simply the product of the 
individual velocity distributions. Moreover, it is neces- 
sary to account for the rate of collisions which depends 
on the relative velocity at the point of impact as well as 
all possibilities of collision by summing up on the loca- 
tions of the impact on the sphere. Pi(z) is then given 
by 

P 1 (z) = cJ J J dndueZv|(v-u).n|/(v) 

/ B (u)<Kz-v.v*), (13) 
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where f(v) and /b(u) are, respectively, the velocity dis- 
tributions of the tagged and bath particles. The integral 
with the subscript S corresponds to the integration over 
the unit sphere with the restriction (u — v).n < where 
n is a unit vector along the axis joining the two centers 
of particles (this imposes that the particles are approach- 
ing each other before colliding). C is the normalization 
constant such that dzP\{z) = 1. 

The postcollisional velocity v* is given by the collision 
rule which is 



1 



rn + M 



[(u- 



(14) 



where M and m are, respectively, the mass of the tagged 
and of the bath particle . a is the normal restitution 
coefficient comprised between and 1. For convenience 0, 
H, 111 , w e introduce a' such that 



1 + a' 



M 



(15) 



A. One dimension 

In one dimension, the integral over angles is replaced 
by counting the right and left collisions. Therefore, 
Eq.JT3J) becomes 



p[ 1 \z) = C / / dudv\u- v\f(v)f B (u)S(z ~ vv*). (16) 



from the case z > 0, where P^ (z) is equal to 



P^(z) 



4C 



(1 + a') 2 

'v(l - a') 



ry/z r+oo 




/ / 


JO Jy/Z 





2z 



1 + a' 



(1 + a')v 



1 )f(v)x 



(20) 



Explicit integration over the velocity can be performed 
and some details of the calculation are given in Appendix 
A. Let us introduce 



M m ( 1 — a 



7 T T \ 1 + a' 



T \ 1 + a! 
2m 1 - a' 
T (1 + a') 2 ' 



(21) 

(22) 
(23) 



the first-collision velocity-correlation distribution p[ X \z) 
then reads for z < 



P^ ] (z)^P^'(0)e 



(24) 



and for z > 



P?\z)=P± ±) (P)e 



+e abz erfc 



a — b 
a + b 
a + b r 



-abz 



1 



(25) 



For granular gases, even if the velocity distributions 
of the tagge d p article and of the bath particle are 
Gaussian [Tol lllj (or close to the Gaussian profile [lj), 
the granular temperatures of these two species are al- 
ways different when a < 1. Let us denote 7 the ratio 
between the bath and the tagged particle temperatures. 
The velocity distribution of the two species read : 



and 



f(v) 



fs(v) 



M 

2nf~ 



-e 2 -t T 



2ttT 



(17) 



(18) 



It is necessary to distinguish the case z < 0, where the 
distribution Pj (z) is given by 



x Jb 



v(l - a') 
1 + a' 



2z 



(1 + a')v 



(19) 



where p[ 1] (0) is given by 



Note that P\ (z) is always asymmetric, contrary to an 
uncorrelated velocity distribution; this results from the 
existence of correlations between pre and post-collisional 

velocities of a particle. Secondly, P x (z) is finite when z 

goes to 0, which is not the case for P£\z). 

To simplify the above expressions, Eq. (|2~I| - (|2~5"f and 
to allow us to discuss the physical results, we now need 
to specify the temperature ratio 7. For inelastic par- 
ticles in a polydisperse granular bath, 7 is in general 
a complicated function of parameters such as the bath 
composition, the heating mechanism and the coefficient 
of restitution. However, three interesting limiting cases 
provide a simple expression of the temperature ratio : (i) 
monodisperse system, for which 7 = 1 (in a Gaussian 
approximation) ; 

(ii) a mixture of granular gases in the limit of infi- 
nite dilution. Martin and Piasecki[l(| showed that the 
velocity distribution of an inelastic tracer in an clastic 
bath remains Gaussian with a granular temperature of 



{a + b){a 2 b 2 -c 2 ) 



2ab\fc 



b 2 -2c 



(26) 
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the tracer given by the relation T e // = "/T^am (equipar- 
tition does not hold); this ratio 7 is given by 



7 



M l + a' 
to 3 — a' 



(27) 



(iii) Thermalized systems of elastic particles (a = 1), 
P^\z) then provides a non trivial information on the 
short-time dynamics even for equilibrium systems. 

In the first case (M = m and 7 = 1). a' = a. This 
gives for the first-collision velocity-correlation distribu- 
tion, for z < 



(1) = M 2 + y/2(l + ^) ^-^^i+^j+i-q), 
1 K ' T 2(1 + a)VTT^ 



(28) 



and for z > 0, 



Jf>(*) 



,1/ 



2(1-°) Ma 
S2 T 



T 2(l + a)Vl + a 2 



(2 - v^(TT^)) e - 



er/ 



y/2- V(l + a 2 ) /M* 



1 + Q- 



r 



+ (2 + v ^TT^)e<^ ( V^))^ 



erfc 



'V2+ v/(l + a 2 ) /Mi 



1 + Q- 



r 



(29) 



The decay of Pj (z) also depends on both the tem- 
perature on the coefficient of restitution a. 

In the second case (granular tracer in a thermalized 
bath), substituting Eq.([27]) in Eqs.([2"T |) -([2"2 ]) . one obtains 
that a = b = jr^p > which gives a simple expression 



for P^\z). For 2 < 



Jf>(*) 



T (1 + a') 3 / 2 



(30) 



and for z > 



T (1 + a') 3 / 2 



erjc 



2mz \ 
—) 

(31) 

In order to show that the local environment influences 
the first-collision velocity-correlation distribution, we re- 
express P^ (z) in terms of the temperature of the tracer 
7T. Substituting Ea.([2?)) in Eas.(|30 j) -([3"T] ) . one obtains 

that for z > P{ >(z) - M 



, e (i+ a >) ,t _ Compar- 



ison with Eq. (f2"8")l shows that P^\z) is sensitive to the 
heating procedure, the temperature of the tracer particle 
being the same in both cases. 

Finally, for elastic hard particles (a = 1), the equipar- 
tition holds, 7=1, and therefore Eqs. (|2"Tj) - (|23[) become 




0.01 



FIG. 1: Log-linear plot of P[ 1] ' (z) versus z (M/T is set to 1) 
(Eqs p8)) - ((29)) ').: Left, from top to bottom, the coefficient of 
restitution a = 0, 0.2, 0.4, 0.6, 0.8, 1. 



a = b = (M + m)/(2VmT) and c = (M 2 - to 2 )/(4toT), 
and P 1 (z) reads for z < 



M 
2T 



M M(M + m ) . 
e 1-m.T ' 

TO 



(32) 



and for z > 



P«(z) 



M 
2T 



M MfM+rnj 

— e 2mT z er/c (M + m) 



2toT 
(33) 



Figure Q] displays P^\z) as a function of z for 
different values of the coefficient of restitution a = 
0, 0.2, 0.4, 0.6, 0.8, 1 for a monodisperse system fEqs. (f28|) - 
([2"9"| ). For elastic hard particles, P^\z) is plotted in 
FigEl as a function of z for different values of the mass 
ratio to/A/ = 1, 1/2, 1/5, 1/10 

Averaged quantities can be deduced from the first- 
collision distributions: The integral of P^ (z) over z < 
corresponds to the fraction of events in which the parti- 
cle velocity after collision is opposite to the precollisional 
velocity. Integrating Eq. (f2~4"|) over all negatives values of 
z leads to 



- , p 1 (1) ( z <o) 



(M 2 + ^7(777 + M)){M 2 - yfy(M - ma)) 
(1 + a)M 2 y / mM + 7TO 2 



with 



M 2 = V(l + a) 2 mM + "f(M - ma) 2 



(34) 



(35) 



For elastic particles (i.e. a = 1), the fraction of colli- 
sions in which the post-collisional velocity has a direction 
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FIG. 2: Log-linear plot of P\ (z) versus z, for elastic hard 
spheres (Eqs.([32])- ([33])). The mass ratios are j% = 1, |, |, ^ 
(on the left side, from top curve to bottom curve). 



FIG. 3: / a), , versus mass ratio m/M for different values 

PI (z<0) ' 

of the coefficient of restitution a: from top to bottom, a = 
0.99,0; 8, 0.5,0. 



opposite to that of the pre-collisional velocity becomes 
very simple because the equipartition property is satis- 
fied (7 = 1): 

V(*<o) = V^W ( 36 ) 

Thus, for a monodisperse system, the probability of hav- 
ing a velocity inverted after a collision is higher than the 
probability of having a velocity whose direction is not 
changed by the collision (in ID). 

For inelastic particles (a < 1), Fig. [3] shows that 
/„(i). „, increases with the mass ratio m/M for all val- 

P{ (z<0) ' 

ues of a. For the sake of simplicity, 7 is assumed equal to 
1. Attention must be paid to the case a = 0: in this case, 
the limit of I p Wr z<0 \ when m/M — > 00 is equal to 1/2 
whereas I p Wr z<0 -\ g° es to 1 when m/M — > 00 if a > 0. 
This discontinuity is clearly apparent in FigJ3] 

B. Two dimensions 

The first-collision velocity-correlation distribution can 
also be obtained analytically in two dimensions and 
above. Indeed the calculation can be performed following 
a method used in the one-dimensional case. Let us note 
that the scalar product of the pre- and post-collisional 
velocities v.v* can be expressed in an orthonormal basis 
associated with the collision as 

v.v* = v n v* n + v t v*, (37) 

where v n and Vt denote the normal and tangential com- 
ponents of the velocity. The post-collisional quantities in 



the rhs of Eg. ([57)) can be eliminated by using Eq. fLij) : 



9 1 — a' , 1 + a' 
v.v* = v 2 t + — 2— < + — 2— v n u n , (38) 



where u n is the normal component of the velocity of the 
bath particle. 

Since the normal and tangent components of the tracer 
particle are independent random variables, the integrals 
can be performed successively over u n and v n as in the 
one-dimensional case, provided that z is shifted to z — vf. 
It is worth noting that the last integration over angular 
variables (n) is trivial, since the integrand does not de- 
pend on n. 

(2) 

For z < 0, the decay of P{ (z) remains exponential 
and is equal to 



P[ 2 \z) = P^(Q) (39) 



For z > 0, P x (z) has two contributions, 
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P^ 1 \o)e ( - ab+c ' >z 
/l + 2(ab + c)^ 



■.erfc 



^ab(z-y) 



erfc 



(a + b)^{z-y) 



V2 



a — b 

r< 

a + b 



-ab(z-y) 



erf 



{a - b)y/(z - y) 
V2 



(a + b)(ab- c) 



1 



Pl {Z<0) M»/a' + P-2cJ 1 + 2iab + c) 



7T 
M 



/ P 1 (2) (2<0) 



m 



M \ M + 2m 



/ P, <1) (z<0) 



(40) 



(2) 

Note that P{ (z) has an exponential decay when the 
scalar product of velocities is negative (i.e. z < 0), with 
the same coefficient {ab + c) than we have obtained in 
one dimension. The integration over vt only changes the 
normalization constant compared to the one-dimensional 
case. Conversely, for positive arguments z, the shape of 

(2) 

P± (z) is more complicated than in one dimension. As 
will be seen in the next section, the exponential decay of 

( 2) 

P{ (z) is universal since this result is the same in any 
dimension and for different heating mechanisms (through 
the 7 dependence). The fraction of events in which the 
scalar product between the pre and post-collisional ve- 
locities is negative is then given by 



(41) 



where a, b and c are given by Eqs. (|2ip 

For instance, for elastic hard particles, the fraction of 
events in which the scalar product is negative after a 
collision is given by 



(42) 



Therefore, for a monodispersc system (m = M), 
-- 1/V6 ~ 0.40824.. in 2D whereas 

1/V2 ~ 0.707.. in ID. In other words, in 
2Z), most collisions do not change the scalar product, 
whereas the converse is observed in ID. 

C. Three dimensions 

In three dimensions, the first-collision velocity- 
correlation distribution can be similarly derived. The 
scalar product between the pre- and post-collisional ve- 
locities can be expressed as 



where u n is the normal component of the velocity of 
the bath particle. P x (z) is the probability distribu- 
tions associated with the sum of two independent ran- 
dom variables. The first one is the ID collision term 
^r-v n u n and the second one is the sum 



'II 



which is a xi-distributed variable, and, as a 



result, an exponentially distributed variable 



, M My 

xi(y) = — e 



27T 



(44) 



(3) 

P| (z) can be expressed as the convolution product 
of the two probability distributions of these variables. 
As the distributions are normalized, the distribution ob- 
tained via the convolution product is normalized. For 

z < 0, 



P["\z) =P«(0) 
=P«(0) 

,(•'5) 



M 



+ DC 



dye 



e (ab+c)z 

l + 2{ab + c)^§ 



-4w e ( ab + c )( z ~y) 



(45) 



and for z > 0, P{ (z) is the sum of several contributions 



p^\z)^pr(o) 



M 



,(ab+c)z er f c 



(a+b)^z~ 
V2 



7 T M, + 2(ab + c) 



M {b~a)e 



erf 



7 T 

(a— b)y/7 
V2 



lT ( a + b)(-±L + 2ab-2c) 



M 



4abe~^* J a 2 + b 2 -2c- ^erf 



yi v /q2 +b 2„ 2c _ g 
72 



7 P (a + b)(-^ + 2ab-2c)(^+2(ab + c)) 



7 T 



(46) 



where a, 6, and c are given by Eq. (f2"Tj) -(f2"3" |) and 7 is the 
temperature ratio. 

As already noted, P^ (z) has the same z dependence 
as in ID and 2D when z < 0. The influence of the di- 
mension is in the amplitude factor which decreases when 
the dimension increases, the temperature and other mi- 
croscopic parameters (masses, coefficient of restitution) 
being kept constant. 

From Eq. (|4"5"|) . the fraction of events with a negative 
velocities scalar product can be exactly obtained; more- 
over, a general formula can be obtained in any dimension: 



v n v n +v t v t + v z v z 



1 



1 



-v n u n + v t +v z 



(43) 



(a + b)(ab - c) 



'n W (*<0) 2abVa^T, 



-2c 



\ + 2{ab- 



(47) 
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IV. INFLUENCE OF THE SONINE 
CORRECTIONS OF THE VELOCITY 
DISTRIBUTION ON P^z) 



For granular gases, whose kinetic properties are well 
described by the Boltzmann equation, the velocity dis- 
tribution is no longer a Gaussian. The deviations from 
Gaussian behavior can be captured by Sonine corrections. 
It is then possible to quantify the influence of these cor- 
rections on the distribution P\{z). (since the definition of 
the latter does not depend on the details of the velocity 
distribution function) . 

The non-linear Boltzmann equation can be solved nu- 
merically by using a Direct Simulation Monte Carlo 
(DSMC) method [la, 0]. We have performed DSMC 
simulations for monodisperse homogeneous systems ex- 
cited through a stochastic thermostat. This allowed us 
to compare the distribution P\{z) obtained by simulation 
with the theoretical prediction calculated with a Sonine 
correction. 

The calculation is similar to that of Eq. (fl"9]) . but in 
the case studied here we no longer consider a tracer par- 
ticle, so that / = Jb- Calculations were performed by 
considering only the first correction. This gives : 



/(v) 



d/2 



M V 

e 



2ttT 



5^- l + a 2 ( a )S- 



/Mv 2 



V 2T 



(48) 

where S 2 is the second Sonine polynomial expressed 
in the appropriate dimension. The value of 02(0;) is 
taken equal to its usual approximation for the stochas- 
tic thermostat 0. [l5j. As the form of the perturbation 
introduced is simply a multiplicative polynomial, the ex- 
pressions obtained for the ID distribution P±(z) are 
analytical. As the result for z > is involved, we will 
only provide that obtained for z < 0, which is 



2(l-Q + N /2v'l + a 2 ) 



if>(z) = 



64V2^(l + a + a 2 + a 3 )" 



x ( Q(a) + a 2 (a) iQ (a) + —Qi(a) + -j^Q 2 (a) 

(49) 

where the <3i(a)'s are simple functions of a given in 
Appendix B. 

Since the 2D case is closer to possible experimen- 
tal systems, we have also included Sonine corrections 

to the velocity distribution for calculating the velocity- 

(2) 

correlation distribution P{ (z), but the lengthy expres- 
sions are not given here. Figure H] displays the analytical 

(21 

result (P{ (z) with Sonine corrections) and the DSMC 
results for two values of the coefficient of restitution : 
a = 0.2,0.5. Even for the more inelastic case, the agree- 
ment between analytical results and DSMC is remark- 
able. 




FIG. 4: DSMC results(symbols) and theoretical predictions 
with the first Sonine correction (lines) for P± 2 \z) versus z = 
rnviv[/T in 2D. Diamonds and circles correspond respectively 
to the simulations for a = 0.2, 0.5, respectively. 



V. SECOND AND HIGHER-COLLISION 
VELOCITY-CORRELATION DISTRIBUTIONS 

The collision statistics can be followed beyond the 
first-collision distribution. Formally, the second-collision 
velocity-correlation distribution P2(z) is given by the re- 
lation 



H - vi).ni 



P2{z) = C2 / / / / GJniG?n 2 duidu2Giv|(ui 
JSi Js 2 J J 

|(u 2 - v*).n|/(v)/ B (u 1 )/s(u 2 ) ( S(z - v.vj), (50) 



where Vi denotes the pre-collisional velocity of the 
tagged particle for the first collision, the velocity af- 
ter the first collision and the post-collisional velocity 
after the second collision. Ui and 112 correspond to the 
velocities of the bath particles for the first and second 
collisions, respectively (Recall that /(v) and /s(v) are 
the velocity distributions of the tagged and bath parti- 
cles, respectively). Obviously, v 2 = v*, and by using 
Eq. (fT4"|) , the collision rule gives for the two collisions 

1 + a' . 

v x =vi-| — [(ui - vij.nijni, (51) 

[(u 2 - v 2 ).n 2 ]n 2 . (52) 



2 

1 + a' 



Finally, C 2 is the normalization constant ensuring that 

/oo 
dzP£\z) = 1. (53) 
-OO 
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In a similar way, it is possible to write down closed equa- 
tions for P„(z). However, whereas tractable expressions 
have been obtained in any dimension for Pi(z), the cal- 
culation increases drastically in complexity for obtaining 
the distribution at the second collision. In the restricted 
case where a' = 1 in one dimension, it is nonetheless pos- 
sible to obtain the exact expression of P 2 (1) (z) (details of 
the calulation are given in Appendix C). Thus, for z < 0, 

^\plabz 



dv- 



f ( hv 

erfc (,71 



erf 



( bz 



\V2v 
(54) 



and for z > 0, 



— y/2abz 




( (ygq - 26)Vi \ 



„ / (-6 + y/a? + fr 2 )y/l \ 



-erfc 



by/a 2 +b 2 z 



-erfc 



2tt(1 + b 2 z) 



V2 

' (b + Va 2 + 6 2 )Vi \ 
V2 ) 



bv 

erfc —= 
V2 



dv- 



erf 



bz 
V2i 



(55) 



where C 2 is determined from Eq. ([55)) and, with the help 
of Eqs.([2T|) and @, a = and b = ^/W. From 



7T """" " V T ' 

Eqs.([54]) and ([55]) . one obtains the small- z expansion for 
P«(z) : 



P 2 (1) W 



V 7 T 



(56) 



Therefore, the second-collision velocity-correlation distri- 
bution shows a divergence at z = 0, reminiscent of the 
divergence of Pd^ (z), when the two velocities are com- 
pletely uncorrelated. However, note that the coefficient 
of the modified Bessel function of the second kind K${z) 

is ab instead of a 2 for P^' (z) . To continue the analysis of 
the velocity-correlation distributions in general, we have 
performed DSMC in various situations, and monitored 
several P„ (z)'s. 

Figure [5J shows Pi 1} (z) as a function of z for n 
i . 2. .10. Note that for n > 3, the distribution practically 
reaches the asymptotic value, namely P oa (z). A simple 
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FIG. 5: nth collision velocity-correlation distributions P^\z) 
versus z for various values of n: n — 1, 2, ..10 (D — 1). 
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FIG. 6: nth collision velocity-correlation distributions P„ ( 
versus z for various values of n: n = 1, 2, .., 10 (D — 2). 



physical interpretation is that after 3 collisions the sys- 
tems loses memory of its initial velocity configuration, 
and the correlation between the initial velocity and the 
velocity after n collisions vanishes when n is larger than 
3. Similar plots are displayed in Figs. [5] and [7] for P„ (z) 

and P^ 3 ' (z), respectively. One notes that the convergence 
to the asymptotic function, P oc (z), becomes slower when 
the space dimension increases. 
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FIG. 7: nth collision velocity-correlation distributions Pi 3 ' (z) 
versus z for various values of n: n = 1, 2, .., 10 (D — 3). 



VI. CONCLUSION 



We have introduced new velocity-correlation distri- 
butions that capture the early stages of the dynamics. 
These non-trivial quantities are efficient probes for inves- 
tigating the environment of a particle in granular gases. 
We have shown that these distributions decay exponen- 
tially when the scalar product of the velocities is negative. 

Other interesting physical situations could be consid- 
ered: free cooling states, in which the velocity distribu- 
tion obeys a scaling form, binary mixtures [la. UtI [HI], 
and not only in the infinite dilution limit considered in 
this paper. 

These distributions are easily accessible in computer 
simulations, and probably, in experiments; indeed, when 
the time resolution is smaller than the mean collision 
time, the probability that two collisions occur during a 
time step is small and the collision history can be moni- 
tored accurately, which permits to build the first-collision 
velocity-correlation distributions . 



APPENDIX A: ID FIRST-COLLISIONS 

We have the following integrals, for z < 0, 
" +0 ° ' aV b 2 z 2 \ a + b Tn 



\ dv(l g ) ex P 

Jo v 



2v 2 



and for z > 0, 



(Al) 



r a ( z u ( aV 



2v 2 



( a -^^ 



(A2) 



and 



f°° . z , 
/ dv(l rjcxp 



t: 1 



2v 2 



(a + b)e abz erfc 



which gives for z < 
P?\z) 



a + b 
~7T 



AC 



(1 + a') 2 ab 



(A3) 



(A4) 



and for z > 



^ [Z) ~ (1 + a') 2 V 2 



a ~ b -abz ( i ( a ~ b 

' erf I — — V 



V2 



a + b abz j- ( a + b 

H —e erfc 



ab 



V2 



(A5) 



The constant C can be obtained by calculating the nor- 
malization condition J dzP\(z) = 1, which gives 



C = 



>Y2 



a 2 b 2 - c 2 



4 x /2ir(a 2 + b 2 - 2c) 



(A6) 



APPENDIX B: SONINE CORRECTION TO THE 
ID CALCULATION 

For z < 0, one finds 
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P^'(z) = e(>+»i 



— T T Q 



du(l - -5) 



ex P ( -jr 



M ({l + a 2 )u 2 



2z< 



(1 + a) 



u 2 (\ + a) 2 , 

Mu 2 (2z/u 2 - (1-a)) 2 ' 



T 



2(1 + a) 2 



(Bl) 
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After integration of Eq. (|Bl[) . one obtains 

P l(0 ) K 1 -^^ 



Q(a) + a 2 (a)Q (a) 

/ Mz M 2 z 2 

Q(a) + a 2 (a) ( Q (a) + — Qi(a) + -j^Q 2 {a) 

(B2) 



where 



Q(a) =16(1 + a) 4 (1 + a 2 )' 

i + a 2 + vWiT 



(B3) 



Qo(a) -3(1 + a) 4 



2 + 6a 2 + 6a 4 + 2a 6 



2Vl + a 2 (l + 6a 2 + Q 4 



(B4) 



Qr(a) = 2^2(1 + a) 2 (l + a 2 ) 3/2 
(37 - 12a + 42a 2 - 20a 3 + 17a 4 ) 
+ 2(1 + a) 2 (1 + a 2 ) (54 - 16a 
+ 84a 2 - 32a 3 + 54a 4 - 16a 5 ), 



(B5) 



where fs(u) denotes the bath velocity distribution and 
f(v) the tagged particle velocity distribution. 

We first integrate on the velocity of the bath particle 
2, namely the velocity of the colliding particle at the 
second collision. We drop the subscript of the velocity of 
tjie bath particle for the collision 1, and P 2 (z) reads 



P 2 (1) (z)=C 2 j dvj du\u-v\\^- U \f B (u)f B (^) /(«). 

(C2) 

Let us introduce the function I(z, v) : 



I(z,v) = / du\u — v\ \z — uv\e 



(C3) 



When z < v 2 , one has 



2w 2 » 2 , 2 2z t 2 , 2 V§k\v\(1 + b 2 z) 



(C4) 



and when z > v 2 



Q 2 {a) = 8 (1 + a 2 ) (3 - 2a + 3a 2 ) 
((7 - 6a + 3a 2 ) (l + a 2 ) + V2\fl + 
(5 - 4a + 5a 2 -2a 3 )). 



(B6) 



APPENDIX C: SECOND-COLLISION 
VELOCITY-CORRELATION DISTRIBUTION 

In one dimension when a' = 1, the velocities of the 
tagged particle and the bath particle are exchanged dur- 
ing the collision. This drastically simplifies the expres- 
sion of the second-collision velocity-collision distribution 
and the calculation becomes tractable. Indeed, if a! = 1, 
P^\z) becomes 

P^\ z ) = ( ^i J dv J dui j du 2 \ui - v\\u 2 - ui|/s(ui)x 
x f B (u 2 )f(v)S(vu 2 - z), (CI) 



e 2„ 



b 2 b 2 
2ttu(1 + b 2 z 



b 3 



erf 



bz 
J2v 



^\v\(l + b 2 z) 
b 3 

Jbv 



(C5) 



Inserting Eq. (|C4[) in Eq. (|C1[) and integrating out the 
first two terms of the integrand leads to Eq. (|5"4"|) . For 
z > 0, by using the property of I(v, z) (Eqs. (|C4|) - JC5|), 

for z > 0, P 2 (z) is expressed as 



p£\ z ) = C 2 I / do - I 
'Vz Jo 



q 2 „ 2 b 2 * 2 

e 2 ~2~J2~ 



I(v,z) 



(C6) 

Integrating out the first terms of the right-hand-side of 
Eq.jC6| leads to Eq.([55]). 
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